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Решена задача о температурном поле нетрадиционного твердого топлива при пиролизе с учетом его структурной неоднороднос- 
ти и неравенства температур газовой и конденсированной фаз твердого топлива во всём диапазоне изменения температур от 
начальной до температур, соответствующих полному завершению процесса пиролиза. Установлено, что перепад температур по 
толщине слоя топлива достигает 600 К при типичных условиях теплообмена, а максимальная разность температур газовой и 
конденсированной фаз слоя твердого топлива составляет 350 К при температуре внешней среды 1000 К. 


Введение 

Теоретический прогноз параметров подготовки 
и сжигания нетрадиционного твердого топлива яв- 
ляется, возможно, одним из важных приложений 
общей теории горения [ 1 ] . Основными элементами 
комплекса физико-химических процессов, пред- 
шествующих воспламенению и горению нетради- 
ционного твердого топлива (НТТ), являются про- 
цессы сушки и пиролиза. В настоящее время эти 
два процесса принято рассматривать в рамках 
приближенных моделей [1, 2]. В основе этих моде- 
лей лежит допущение о возможности осреднения 
по толщине слоя топлива значений всех искомых 
функций и характеристик и решения задач сушки и 
пиролиза с использованием систем обыкновенных 
дифференциальных нелинейных уравнений [2] или 
даже более простых математических моделей. 

В то же время известно, что в аналогичных зада- 
чах теории тепловой защиты используют модели, 
учитывающие, не только неоднородность основ- 
ной искомой функции температуры по простран- 
ственной координате [3], но также и неравенство 
температур газообразных и конденсированных 
продуктов реакции термического разложения теп- 
лозащитных материалов [4, 5]. Поэтому целесооб- 
разным является анализ масштабов возможной 
разности температур газовой и конденсированной 
фаз слоя НТТ в режимах сушки и пиролиза. 

Цель данной работы - моделирование процесса 
пиролиза нетрадиционного твердого топлива, а так- 
же оценка масштабов возможной пространственной 
неоднородности температуры слоя НТТ как по ко- 
ординате, ориентированной по направлению нор- 
мали к нагреваемой поверхности, так и по коорди- 
нате, параллельной нагреваемой поверхности. 

Постановка задачи 

Рассматривается задача о температурном поле 
слоя НТТ (опилок или отходов деревообработки) в 
форме прямоугольной пластины бесконечной дли- 
ны, нагреваемой по верхней грани за счет конвек- 
тивного и радиационного теплообмена с внешней 
средой. При нагреве НТТ до температуры начала пи- 
ролиза Т ш = 473 К начинается термическое разложе- 
ние топлива, протекающее [1] с образованием кон- 
денсированных и газообразных продуктов реакции. 


Принята следующая геометрическая модель ис- 
ходного топлива и образующегося в процессе пи- 
ролиза кокса. Нетрадиционное твердое топливо в 
исходном состоянии представляет собой неодно- 
родную по своему составу пористую среду. Эта сре- 
да обладает определенной начальной пористостью. 
Нижняя граница области решения задачи является 
границей слоя. Правая и левая границы области ре- 
шения совпадают с осями симметрии газовой и 
конденсированной фаз НТТ (рис. 1, о). 



Рис. 1 . Схема области решения задачи: 1) исходное топливо; 
2) НТТ в состоянии термического разложения ; 3) кокс 
НТТ (процесс пиролиза завершен); 4) стенка; 5) газо- 
вая фаза слоя НТТ (исходный воздух и газообразные 
продукты разложения); (а) — начальное состояние; 
(Ь) — промежуточное состояние 
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Эволюция структуры НТТ при пиролизе предс- 
тавлена схематично на рис. 1 , Ъ. По мере продвиже- 
ния зоны пиролиза вглубь топлива образуется кокс 
[1]. При этом изменяется исходная пористая струк- 
тура за счет увеличения размеров начальных пор 
естественного происхождения и изменения их 
формы в результате пиролиза. Для структуры слоя 
НТТ как в начальном состоянии, так и после завер- 
шения процесса пиролиза, можно выделить два ха- 
рактерных размера: толщину конденсированного 
вещества Д к и расстояние между стенками пор Д г . 
Сумма Д к + Д г = сі = соп8І, характеризует структуру 
и пористость НТТ. Для значений Д к была принята 
следующая аппроксимация 

Л к = Лко[1-(1-ад при Т ш <Т<Т т , X = Р ■ 

Ро Рк 

Здесь Т ш и Т ш - температуры начала и конца ре- 
акции пиролиза соответственно; К- коксовое чис- 
ло; х ~ степень разложения; Д ко - начальное значе- 
ние Д к (при Т < Т ш )\ р 0 , р, р к - начальная, текущая 
и конечная плотности НТТ, соответственно. 

Принятая геометрическая интерпретация слоя 
НТТ достаточно полно, по сравнению с моделями, 
использующими осреднение всех искомых функ- 
ций и параметров по толщине слоя, описывает тем- 
пературное поле топлива с учетом его реальной 
структуры. Характерный поперечный размер й, 
описывающий как исходную неоднородность НТТ, 
так и изменение структуры в процессе пиролиза по 
данным [1,2] изменяется от 1 до 2 мм. В этом диа- 
пазоне изменялось значение сі и при проведении 
численных исследований. 

Задача решалась при следующих основных ог- 
раничениях и допущениях: 

1. Процесс термического разложения топлива мо- 
жет быть описан одной эндотермической реак- 
цией [1]. 

2. При нагреве и термическом разложении НТТ 
не происходит "вспучивания" [6] или усадки, а 
также отсутствуют механические деформации 
материала, обусловленные внешними воздей- 
ствиями или массовыми силами. 

3. Лучистый теплообмен в слое кокса по своим 
масштабам значительно уступает теплообмену 
за счет конвекции и кондукции и может не учи- 
тываться. 

4. Задача рассматривается для двухслойной плас- 
тины конечных размеров "НТТ-стенка". Тол- 
щина стенки при этом выбирается близкой к 
толщине слоя НТТ, и на внешней ("холодной") 
границе стенки выставляется условие теплоизо- 
ляции. 

5. Топливо в исходном состоянии представляет 
собой пористую структуру. В процессе пироли- 
за по мере продвижения фронта начала реакции 
вглубь материала, конденсированные продукты 
реакции стягиваются в монолитную углеродис- 
тую пластинку (стенку поры) для каждого эле- 
ментарного характерного размера сі. 


6. Не учитываются возможные газофазные хими- 
ческие реакции между продуктами пиролиза, а 
также возможное осаждение пироуглерода на 
стенках пор. 

7. Масштабы конвективного теплопереноса в по- 
ристой структуре топлива определяются массо- 
вым расходом газообразных продуктов пироли- 
за гіі. 

8. Твердое топливо не содержит влаги (процесс 
сушки завершен). 

9. Теплообмен на границе "НТТ-стенка" прини- 
мается идеальным. 

10. Концентрация окислителя (воздуха) на нагре- 
ваемой границе слоя НТТ недостаточна для его 
воспламенения при температурах до 1000 К рав- 
ной температуре конца пиролиза Т ш . 


Математическая модель 

С учетом принятых допущений процесс тепло- 
массопереноса в неоднородной пористой структу- 
ре НТТ и стенке (рис. 1, б) может быть описан с ис- 
пользованием системы уравнений, краевых и на- 
чальных условий, вытекающих из общей математи- 
ческой модели газификации твердого топлива [1]. 
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д з ІТа = Я 5 гТа, Т ъ =Г 5 ,ж = Д„ 0<у<у т ; 
дх дх 
дТ 2 дТ 5 

= 4; Т = * = д , > Л, < X ^ Л, ; 

дх дх 

Т 1 ^- = Я 5 ^-,Т=Г 5 ,Л- = Д, 0! ѵ ш < я < 5 ; 

I = О, Т, =Т 2 =Т Ъ =Т Л =Т 5 =Т 0 , у нп = О, у т = 0, Д„ = Д„ 0 , Я = 0. 

Здесь: Р - давление; К - универсальная газовая 
постоянная; ср - пористость; М -молекулярная мас- 
са; р- плотность, с -удельная теплоемкость, (- вре- 
мя; х, у - прямоугольные координаты, Я - коэффи- 
циент теплопроводности, т - массовый расход газо- 
образных продуктов пиролиза исходного материала, 
\Ѵ и 0 - скорость и тепловой эффект реакции пиро- 
лиза; а - постоянная Стефана-Больцмана; е щ - при- 
веденная степень черноты; а - коэффициент тепло- 
обмена; 5 - толщина слоя НТТ; к - толщина стенки; 
с! - ширина рассматриваемой области. 

Индексы: 0 - начальное значение; в - внешняя 
среда; нп - начало пиролиза; кп - конец пиролиза; 
ст - стенка; 1 - исходный НТТ; 2 - топливо в сос- 
тоянии термического разложения; 3 - кокс (про- 
цесс пиролиза завершен); 4 - стенка; 5 - газовая 
фаза топлива (исходный воздух и газообразные 
продукты разложения НТТ). 

Метод решения 

Решение сформулированной задачи проведено 
методом конечных разностей [7]. Разностные ана- 
логи дифференциальных уравнений решены ло- 
кально-одномерным методом [7]. Решение разно- 
стных одномерных уравнений проведено методом 
итераций с применением метода прогонки на каж- 
дой итерации [7]. Использовалась неравномерная 
разностная сетка. Численный анализ выполнен при 
следующих значениях параметров и постоянных: 
Т ъ = 800... 1200 К; а = 30 Вт/м 2 -К; А 3 = 1,0 Вт/м-К; 
с 3 = 1200 Дж/кг-К; с 5 = 1200 Дж/кг-К; р 3 = 1 600 кг/м 3 ; 
<5 = 1...4Т0 2 м; к = (10...30) 10 2 м; щ = 0,1...0,3; 
Я; = 0,1 Вт/м-К; С; = 1500 Дж/кг-К; р, = 650 кг/м 3 ; 
Я 4 = 0,6 Вт/м-К; с 4 = 780 Дж/кг-К; р 4 = 1500 кг/м 3 . 

Значения Я 5 газообразных продуктов пиролиза 
при температурах от 400 до 1100 К рассчитывались 
с использованием табличных данных [8] для газов, 
близких по своему составу к продуктам пиролиза 
древесины. Кинетические постоянные пиролиза 
принимались соответствующими энергии актива- 
ции и предэкспоненте для пиролиза древесины [9]. 

Наличие подвижных границ по координатам х и 
у обуславливало необходимость использования не- 
регулярной и неравномерной разностной сетки, 
параметры которой выбирались из условий обеспе- 
чения сходимости и устойчивости разностного ре- 
шения сформулированной задачи. Перестройка 
сетки проводилась на каждом временном слое. 
Вблизи границ раздела газовой и конденсирован- 
ной фаз проводилось сгущение сетки (увеличения 
числа узлов по пространственным координатам на 
единицу длины по сравнению с областями газовой 
и твердой фаз). Значения температуры и других ис- 


комых функций в узлах новой сетки определялось 
интерполяцией по значениям функций в узлах ста- 
рой сетки. Общее число узлов при этом составляло 
от 50000 до 200000. 

Численное исследование сходимости проводи- 
лось на множестве сеток с уменьшением значений 
шагов по пространственным координатам. Прово- 
дилось также численное исследование сходимости 
итерационного процесса. Были получены оценки 
равномерной нормы разности между значениями 
температуры, соответствующими двум последова- 
тельным итерациям. Порядок невязок для типич- 
ных условий решения задачи не превышал 0,1%. 


Результаты и обсуждение 


На рис. 2, 3 приведены типичные результаты 
расчетов температурных полей НТТ в процессе пи- 
ролиза в условиях теплообмена, соответствующего 
типичному режиму на внешней (нагреваемой) гра- 
нице слоя материала. 




Рис. 3. Распределения температуры по координате у в двух 
сечениях по х в различные моменты времени: 1, 3, 5) 
конденсированная фаза НТТ; 2, 4, 6) газовая фаза 
НТТ; 1, 2)1 = 20 с; 3, 4) I = 90 с; 5, 6) { = 300 с; 
(Т, = 1000 К; ср 50 = 0, 1; 8-0,04 м; Я 3 = 1,0 Вт/м-К) 

Распределения температуры по координате у в 
двух сечениях по х, проходящих через оси симметрии 
газовой и конденсированной фаз одиночной поры, 
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показывают, что разность температур газообразных и 
конденсированных продуктов пиролиза НТТ дости- 
гает после 20 с начала процесса почти 350 К при тем- 
пературе поверхности около 1000 К. С ростом време- 
ни разность температур несколько уменьшается, но 
даже при /=300 с, когда фронт пиролиза достигает 
границы стенки, составляет около 150 К. При этом 
более высокой является температура газовой фазы. 

Полученный результат может быть объяснен бо- 
лее высоким (почти в 10 раз) коэффициентом тем- 
пературопроводности газа по сравнению с материа- 
лом твердой фазы продуктов пиролиза НТТ и отно- 
сительно низкой ролью конвективного теплообмена 
газообразных продуктов пиролиза со стенками пор. 

Следует особо отметить, что перепад температур 
по толщине слоя НТТ для типичных режимов [1] 
составляет почти 600 К для интервала времени от 
нуля до 300 с. Такой большой градиент температур в 
интервале времени, представляющем наибольший 
интерес, возможно должен соответствовать интен- 
сификации процесса пиролиза и, соответственно, 
горения НТТ, по сравнению с моделями, использу- 
ющими осредненные по толщине слоя характерис- 
тики (в том числе и температуру). Разный вид кри- 
вых Т(у) на рис. 3 для газовой и конденсированной 
фаз обусловлен различным характером изменения 
теплофизических характеристик этих фаз НТТ в 
состоянии пиролиза. Увеличение значения Я 3 с рос- 
том температуры приводит (при переходе материа- 
ла частицы древесины в состояние углеродистой 
стенки поры кокса) к перегибу на кривой Т(у) для 
конденсированной фазы. Изменения Я 5 с ростом Т 
менее значительны и не играют такой заметной ро- 
ли в формировании температурного поля газооб- 
разных продуктов пиролиза топлива. 

Полученные результаты численного анализа 
позволяют сделать вывод о возможной целесооб- 
разности использования в дальнейшем при анали- 
зе процессов тепломассопереноса, протекающих в 
слое нетрадиционного твердого топлива при пиро- 
лизе, двумерные модели. 

Проведен анализ влияния ряда внешних факто- 
ров и характеристик топлива на масштабы выде- 
ленных особенностей в температурных полях слоя 
НТТ при пиролизе. 

Численный анализ проведен также и для различ- 
ных толщин слоя НТТ. Результаты расчетов (рис. 4) 
показывают, что увеличение <5 приводит при прочих 
адекватных условиях к незначительному росту тем- 
пературной неравновесности газовой и конденсиро- 
ванной фаз, причем этот рост можно отметить толь- 
ко в области относительно низких температур 
( Т <Т Ш ), и он обусловлен, в основном, влиянием 
стенки. Так распределения значений АТ=Т 5 -Т 3 по у 
для четырех различных толщин слоя топлива (рис. 4, 
Ъ ) показывают роль теплоотвода в стенку. 

Кривая 1 рис. 4, Ь соответствующая минимальной 
толщине слоя топлива 10 мм, показывает, что тепло- 
отвод в стенку приводит к такому снижению темпе- 
ратуры газа в порах, что вблизи границы "НТТ-стен- 


ка" Т 5 значительно меньше 7). Увеличение толщины 
слоя НТТ приводит при прочих адекватных условиях 
к появлению локальных минимумов и максимумов 
на кривых Д Т(у). 




Рис. 4. Распределения температуры (а) и разности темпера- 
тур АТ=Т—Т К (Ь) по координате у при различных тол- 
щинах слоя НТТ: (а) 1, 3) конденсированная фаза НТТ; 
2, 4) газовая фаза НТТ; 1,2)8= 0,01 м; 3,4)8 = 0,03 м 
(1 = 180 с; <р Б 0 = 0,1; Т, = 1000 К; Я 3 = 1,0 Вт/м- К); (Ь) 
1)8 = 0,01 м; 2)8 = 0,015 м;3) 8 = 0,02 м; 4)8 = 0,04 м 
(1 = 240 с; <р 50 = 0, 1; Т в = 1000 К; Я 3 = 1,0 Вт/м- К 

Наличие мелкомасштабных локальных миниму- 
мов и максимумов на каждой из трех кривых (рис. 4, Ъ) 
обусловлено как изменением теплофизических ха- 
рактеристик конденсированной фазы при переходе 
материала древесины в науглероженное состояние, 
так и реализацией в определенном диапазоне темпе- 
ратур (от Т ш до Т ш ) процесса пиролиза, протекающе- 
го с поглощением тепла в конденсированной фазе и 
выделением газообразных продуктов. Так, при сме- 
щении "холодной" (ненагреваемой) границы (где 
температуры 7) и Г 5 в течение достаточно длительно- 
го промежутка времени равны) за счет передачи теп- 
ла от нагреваемой границы газ, заполняющий исход- 
ные поры НТТ, прогревается (вследствие более высо- 
кой температуропроводности) быстрее древесины. 
Теплообмен между газом, заполняющим поры, и час- 
тицами древесины в порах НТТ до начала процесса 
пиролиза можно обоснованно считать незначитель- 
ным. По мере смещения к нагреваемой поверхности 
эта разность АТ= Т 5 -Т } растет по указанным выше 
причинам. При достижении фронта начала пиролиза 
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с Т=Т т начинается процесс пиролиза, протекающий 
с поглощением тепла и выделением газообразных 
продуктов. Последнее приводит к некоторому сни- 
жению значений АТ (отмеченному на всех кривых 
рис. 4). К этим же результатам приводит и рост Я 3 , 
обусловленный пиролизом - растет коэффициент 
температуропроводности а 3 . По мере завершения пи- 
ролиза снижение АТ с уменьшением расстояния до 
нагреваемой поверхности прекращается и разность 
температур газовой и конденсированной фаз кокса 
НТТ опять растет, достигая максимальных значений 
вблизи нагреваемой поверхности. В связи с этим сле- 
дует отметить, что учет радиационного теплоперено- 
са в порах, возможно, привел бы к изменению вели- 
чины АТ Но решение аналогичной задачи с учетом 
радиационного теплообмена в пористом слое кокса 
топлива сопряжено как с существенным усложнени- 
ем математической модели, так и с использованием 
ряда эмпирических характеристик радиационных 
свойств кокса. Последние пока известны с недоста- 
точной для детального теоретического анализа на- 
дежностью, поэтому такой анализ пока невозможен. 



Рис. 5. Распределения разности темпера тур АТ по координа - 
те у при различных температурах внешней среды Ту 
1) 800 К; 2) 1000 К; 3) 1200 К; (1=45 с; <р 50 = 0,1 ; 
5=0,04 м; й.з=1, 0 Вт/мК) 



Рис. 6. Распределение разности температур АТ по координа - 
те у при различных коэффициентах теплопроводнос- 
ти А 3 : 1) 0,4; 2) 1,0; 3) 1,5 Вт/мК; (1=60 с; <р, 0 =О,1; 
Т,=1000 К; 5=0,04 м) 

На рис. 5 приведены распределения разности 
температур А Т газовой и конденсированной фаз по у 
в тех же сечениях по х, что и на рис. 3, полученные 
при различных значениях температуры внешней 


среды на нагреваемой границе слоя НТТ. Видно, что, 
как и можно было ожидать, перепад А Г по толщине 
слоя увеличивается с ростом Т в . При этом распреде- 
ления АТ по у имеют немонотонный характер. 

На рис. 6 приведены зависимости АТ(у) для 
трех разных значений коэффициентов теплопро- 
водности твердой фазы кокса НТТ Я 3 . Хорошо вид- 
но, что увеличение Л 3 приводит к деформации кри- 
вых и появлению локальных максимумов и мини- 
мумов разности АТ. Увеличение Я 3 вызывает сни- 
жение АТ в области пиролиза. 

Увеличение же толщины стенки, прилегающей 
к НТТ, практически незначительно сказывается на 
температурном поле НТТ. Возможно, данный ре- 
зультат обусловлен относительно низкой темпера- 
туропроводностью материала стенки. 

В таблице приведены результаты численного 
анализа влияния начальной пористости НТТ на 
температурное поле. Хорошо видно, что увеличе- 
ние начальной пористости в 2,5 раза приводит к из- 
менению температуры газовой фазы менее чем на 
5 %. Данный результат косвенным образом также 
подтверждает незначительное влияние конвектив- 
ного теплообмена в порах на формирование темпе- 
ратурного поля слоя кокса топлива. 

Таблица. Распределение температуры Г 5 (К) по у при 1=120 с 
при разных значениях начальной пористости слоя НТТ 


<йю 

у, ММ 

1 

3 

5 

7 

10 

12 

0,1 

745 

676 

603 

521 

349 

310 

0,25 

754 

684 

613 

533 

361 

317 


Выделенные в результате выполненного теоре- 
тического анализа закономерности процесса пиро- 
лиза древесины при нагреве до высоких температур 
создают предпосылки для использования более 
сложных моделей пиролиза НТТ [1]. Так, напри- 
мер, существенные градиенты температур по тол- 
щине рассматриваемого слоя должны приводить к 
интенсификации процесса распространения тепла 
вглубь топлива и, соответственно, к росту скорости 
пиролиза последнего. 

Заключение 

На основании выполненного численного ана- 
лиза можно сделать вывод о существенной неодно- 
родности температурного поля слоя нетрадицион- 
ного твердого топлива при пиролизе в условиях 
теплообмена на внешней нагреваемой границе, со- 
ответствующих условиям работы систем подготов- 
ки. При этом масштабы отклонения температур га- 
зообразных и конденсированных продуктов пиро- 
лиза древесных отходов достаточно существенны - 
достигают 300 К в различные моменты времени. 
Полученные результаты дают также основания для 
вывода о возможной целесообразности учета тем- 
пературной неоднородности слоя НТТ при реше- 
нии задач о зажигании таких топлив. 


94 





Технические науки 


СПИСОК ЛИТЕРАТУРЫ 

1. Основы практической теории горения / В. В. Померанцев, 
М.К. Арефьев, Д.Б. Ахметов и др. — Л.: Энергия, 1973. —227 с. 

2. Канторович Б.В. Основы теории горения и газификации твер- 
дого топлива. — М.: Изд-во АН СССР, 1958. — 194 с. 

3. Полежаев Ю.В., Юревич Ф.Б. Тепловая защита. — М.: Энер- 
гия, 1976. —349 с. 

4. Ким Л. В. Определение коэффициента теплообмена в порис- 
тых средах // Инженерно-физический журнал. — 1993. — Т. 65. 
— № 6. —С. 663—667. 

5. Алифанов А.М., Трянин А.П., Ложкин А.Л. Эксперименталь- 
ное исследование метода определения коэффициента внутрен- 
него теплообмена в пористом теле из решения обратной зада- 


чи // Инженерно-физический журнал. — 1987. — Т. 52. — № 3. 
— С. 461—469. 

6. Страхов В. Л., Гаращенко А. И., Рудзинский В.П. Расчет неста- 
ционарного прогрева многослойных огнезащитных конструк- 
ций // Вопросы оборонной техники. — 1944. — Сер. 15. — 
Вып. 1(109-110). -С. 30-36. 

7. Самарский А.А. Теория разностных схем. — М.: Наука, 1983. — 
354 с. 

8. Корчунов Ю.Н., Тюльпанов Р.С. Исследование скорости тер- 
мического разложения древесины и торфа // Инженерно-фи- 
зический журнал. — 1960. — № 7. — С. 102—105. 

9. Варгафтик Н.Б. Справочник по теплофизическим свойствам 
газов. — М.: Физматгиз, 1963. —142 с. 


ѵпк чдт е, пн 




